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ABSTRACT 

We present optical and near-infrared colour maps of the central regions of bulges of SO 
and spiral galaxies obtained with WFPC2 and NICMOS on the HST. By combined 
use of HST and ground-based data the colour information spans from a few tens of 
pc to a few kpc. In almost all galaxies the colour profiles in the central 100-200 pc 
become more rapidly redder. We attribute the high central colour indices to a central 
concentration of dust. We infer an average extinction at the center of Ay = 0.6 - 
1.0 mag. Several objects show central dust rings or disks at sub-kpc scales similar to 
those found by others in giant ellipticals. For galactic bulges of types SO to Sb, the 
tightness of the B — I vs I — H relation suggests that the age spread among bulges 
of early type galaxies is small, of at most 2 Gyr. Colours at 1 R e ff, where we expect 
extinction to be negligible are similar to than those of elliptical galaxies in the Coma 
cluster, suggesting that these bulges formed at the same time as the bright galaxies in 
Coma. Furthermore the galaxy ages are found to be independent of their environment. 
Since it is likely that Coma was formed at redshift z > 3, our bulges, which are in 
groups and in the field, also must have been formed at this epoch. Bulges of early- 
type spirals cannot be formed by secular evolution of bars at recent epochs, since such 
bulges would be much younger. There are three galaxies of type Sbc and later, their 
bulges are younger and could perhaps arise from secular evolution of transient bars. 
Our results are in good agreement with semi- analytic predictions (Baugh et al. 1996, 
1998, Kauffmann 1996), who also predict that bulges, in clusters and in the field, are 
as old as giant ellipticals in clusters. 
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1 INTRODUCTION 

We do not have a clear picture of the formation mecha- 
nism for the central bulges of spiral galaxies. Their structural 
and dynamical properties have long suggested the view that 
bulges are like small elliptical components residing in the 
center of a large disk. Indeed, bulges obey the D n - a rela- 
tion and fall on the Fundamental Plane of elliptical galax- 
ies (Dressier et al. 1987, Bender et al. 1992). Their surface 
brightness profiles fall off as fi(r) oc r '" (Andredakis et 
al. 1995), with n = 4 for early- type galaxies falling to n = 
1 for Sc and later types. Dynamically, bulges are consistent 



with oblate, isotropic models, like low-luminosity ellipticals 
(Kormendy & Illingworth 1982, Davies et al. (1983)). In 
a series of papers Peletier and Balcells have studied optical 
and near-infrared (NIR) colour distributions of a complete 
sample of spiral bulges, using ground-based data (Balcells & 
Peletier 1994, Peletier & Balcells 1996, 1997) with the goal 
of extending the comparison of bulges and ellipticals to their 
stellar populations. They find that, once dust is accounted 
for, the colours of bulges are never redder than those of el- 
liptical galaxies of the same luminosity. Bulge colours are 
very similar to those of the inner disk, the differences being 
much smaller than colour differences from galaxy to galaxy. 
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Population models then suggest that the inner disk (at 2 
scale lengths) must have formed at the same time as, or at 
most 3 Gyr after, the bulge, and that bulge metallicities are 
lower than those of giant ellipticals. 

HST allows us to investigate the centres of bulges with a 
ten fold increase in angular resolution, yet HST data on the 
colours of galaxy centres at scales of tens of parsecs is scarce. 
Colours contain useful information on the stellar populations 
of the galaxy centres, and allow us to estimate the level of 
internal extinction. HST data have shown that dust patches 
are very common in galaxy nuclei, even in ellipticals (van 
Dokkum & Franx 1995) , hence we expect a strong signature 
of dust in spiral bulges. NICMOS allows an unprecedented 
look at the inner structure of bulges in the H band where 
the extinction is approximately a factor 6 lower than in V. 
NICMOS can produce more reliable measurements of sur- 
face brightness profiles and isophotal shapes. 

Here we show the results of the first colour study of the 
centers of galactic bulges based on HST/NICMOS data. We 
observed in the NICMOS F160W band (tf-band), F814W 
(7) and in F450W, a wider version of the B-band. This com- 
bination of passbands is very suitable to study stellar pop- 
ulations in galaxies with limited amounts of extinction and 
recent star formation, like centres of early-type spirals. In 
those galaxies the observed colours are determined by the 
underlying old stellar population, perturbed by dust extinc- 
tion and the light of young, recently formed stars. Extinction 
by dust affects all optical and near-infrared colours, while re- 
cent star formation can only be seen in blue passbands (see 
e.g. Knapen et al. 1995). In that paper it is shown that a 
red colour, like I — H will primarily show the extinction due 
to dust, as well as old stellar population gradients, while a 
colour which includes a blue band like B or U, in this pa- 
per B — I, is specially sensitive to recent star formation. In 
the absence of extinction, B, I and H make it possible, in 
principle, to measure the age and metallicity of the stellar 
population (e.g. Aaronson et al. 1978; Peletier et al. 1990; 
Bothun & Gregg 1990). We infer large amounts of centrally 
concentrated dust within the inner 100 pc. Outside this re- 
gion, the colour-colour distribution is quite tight, allowing 
us to place limits on the age dispersion of galactic bulges. 
With the HST it is easier to reach a high photometric accu- 
racy, because there is no atmospheric extinction to correct 
for, the instrumental PSF is stable, and the instrument is 
very well characterized, as compared to ground-based instru- 
ments. This means that we can expect accuracies of 0.02- 
0.03 mag or better in optical and near-infrared colours. Such 
high quality measurements are vital if we are to separate the 
effects of age and metallicity using colours, even though the 
large wavelength baseline (B-I-H) minimises the effect of 
photometric errors. 

Currently there are three main scenarios for the for- 
mation of bulges: bulges form before, contemporaneously 
with, or after disks (see the review by Wyse et al. 1997, 
and Bouwens et al. 1999). In the first scenario (Eggen et al. 
1962, Larson 1975, Carlberg 1984) the formation of bulges is 
mainly described by collapse of a primordial gas cloud into 
clumps, which then merge together. The disk only forms 
after the last massive merger via gas infall. In a second sce- 
nario, infall of a gas-rich dwarf galaxy onto a disk produces 
star formation more or less at the same time in bulge and 
disk (Pfenniger 1992). In the third scenario of secular evolu- 



tion of disks (e.g Combes et al. 1990, Pfenniger & Norman 
1990, Norman et al. 1996), the bulge is formed by dynamical 
instabilities of the disk, in which a bar is formed, which then 
forms a massive central concentration by allowing material 
to stream towards the central regions. The central mass con- 
centration itself then tidally disrupts the bar and forms the 
bulge. When material is continuing to fall onto the disk 
sometime later a new bar can be formed, and the whole 
process will start again and induce more star formation in 
the bulge. To distinguish between these models one needs to 
measure the ages of both bulges and disks. Peletier & Bal- 
cells (1996) and Terndrup et al. 1994, using ground-based 
data, found that the ages of bulges and inner disks are very 
similar. Here in this paper we measure the ages of the bulges 
themselves. 

The paper is organized as follows: Section 2 gives details 
of the observations and data reduction. Section 3a discusses 
the dust content of the galactic bulges. Section 3b describes 
the analysis of the stellar populations, and in Section 4 some 
implications for our understanding of galactic bulges are dis- 
cussed. 



2 SAMPLE AND OBSERVATIONS 

Twenty galaxies were observed with HST in Cycle 7 (Sum- 
mer 1997) with WFPC2 (F450W and F814W) and NICMOS 
(F160W). They all are part of the original sample of Balcells 
& Peletier (1994), which is a complete, B-magnitude- limited 
sample of early-type spirals (type SO-Sbc) with inclinations 
larger than 50°, and for which one side of the minor axis 
colour profile is approximately featureless as seen from the 
ground. The subsample considered here was chosen to in- 
clude galaxies of types SO-Sbc, excluding galaxies that are 
exactly edge-on, and is biased towards the nearest objects. 
The galaxies observed are given in Tables 1 and 2. Their 
basic parameters: luminosities, effective radii and ground- 
based colours, are given in Andredakis et al. (1995) and 
Peletier & Balcells (1997). The mean absolute i?-band mag- 
nitude of the sample is Mr = -22.50 mag, with an RMS 
dispersion of 1 mag (Ho = 70 km s _1 Mpc -1 ). 

The standard HST-pipeline data reduction was used for 
the WFPC2 optical data. For NGC 7457 the /-band ob- 
servations were taken from the HST Archive. For the H- 
band two MULTIACCUM NIC2-exposures of in total 256 
s were taken, offset from each other by 1", to enable us to 
remove bad pixels. The standard STSDAS CALNICA re- 
duction package was used, after which a small fraction of 
the flatfield was subtracted from the reduced images, to ac- 
count for a non-zero pedestal level in the raw image. This 
fraction was determined by requiring that the final image 
was smooth, and amounts to about 0.3-0.5 ADU s , cor- 
responding to about H—17 mag (arcsec) -2 . Applying this 
step is only important for the lowest flux levels, but the 
images improve considerably in quality (they look like the 
ground-based images). The final sky background level was 
determined by comparing the image with the ground-based 
image in the same band. To compare the H-b&nd data with 
the ground-based data in K, we used the following transfor- 
mation, derived from the new GISSEL96 models of Bruzual 
& Chariot (see Leitherer et al. 1996): 



H-K = 0.111(7 - K) 



0.0339 
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Table 1. Basic Parameters 



THE B - I VS. / 
DIAGRAM 



H COLOUR-COLOUR 



NGC 


T 


R c ff 


-£* -*ccn 


B - ha 


-* "ccn 


I — Hcti 


5326 


1 


2.16 


2.510 


2.089 


2.117 


1.796 


5389 





4.67 


2.992 


2.082 


2.208 


1.823 


5422 


-2 


4.51 


2.581 


2.092 


2.170 


1.809 


5443 


3 


4.70 


2.894 


1.941 


2.408 


1.554 


5475 





2.58 


2.342 


2.065 


2.014 


1.774 


5577 


4 


4.14 


1.910 


1.862 


1.792 


1.665 


5587 





2.13 


2.600 


1.993 


2.086 


1.840 


5689 





3.72 


2.355 


2.058 


2.231 


1.816 


5707 


2 


3.40 


2.720 


2.069 


2.144 


1.722 


5719 


2 


3.88 


3.494 


2.252 


2.758 


1.972 


5746 


3 


14.6 


4.117 


2.040 


2.910 


1.749 


5838 


-3 


7.50 


2.600 


2.091 


2.184 


1.767 


5854 


-1 


6.03 


1.872 


1.840 


1.821 


1.659 


5879 


1 


2.29 


2.653 


1.932 


2.353 


1.704 


5965 


3 


8.23 


2.946 


2.051 


2.406 


1.755 


6010 





2.43 


2.543 


2.086 


2.105 


1.843 


6504 


2 


4.98 


2.452 


2.064 


2.062 


1.764 


7331 


3 


5.0 


2.356 


2.047 


1.978 


1.766 


7457 


-3 


8.55 


1.951 


1.920 


1.823 


1.620 


7537 


4 


1.48 


2.928 


1.948 


2.560 


1.969 


4472 


- 


- 


- 


2.357 


- 


1.947 



Notes to Table 1: The morphological types in column (2) are 
from the RC3 (dc Vaucouleurs ct al. 1991). The effective radii of 
the bulges in column (3) in arcseconds are from Andrcdakis et al. 
(1995), except for NGC 7331, where the value is an estimate based 
on Prada ct al. (1996). Distances in kpc are given in Andrcdakis 
ct al. (1995). Also given arc the colours of NGC 4472 at 5" (sec 
text) . 



3.1 Dust in Nuclei of Galactic Bulges 

In Fig. 1 we show intensity-maps of the inner regions of the 
galaxies, together with B — I and I — H colour maps with the 
same scale and orientation, superimposed on //-band con- 
tour maps. The maps indicate the position of dust lanes. It is 
clear from the figure that many, maybe all, galaxies have red 
nuclei. To show this more clearly we show all colour profiles 
obtained in wedges along the dust free side of the minor axis 
in Fig. 2. This is not the first time that red nuclei have been 
found in spiral galaxies. For example, in our ground-based 
data (Peletier & Balcells 1997) we show several galaxies with 
red nuclei. For that sample however the colour profiles on the 
least dusty side are generally featureless at radii larger than 
1", showing the same logarithmic colour gradients through- 
out the bulge. Inside 1" no information is available about the 
colours from ground-based data. The fact that the PSF of 
the HST data is well known and stable enables us to correct 
the colour-profiles and maps for most of the instrumental 
PSF effects. In this way one can measure colours down to 
radii as small as ~ 0.10" (the diffraction limit is 0.15"). At 
this resolution the red nuclei are easy to see, and extended 
in most of our sample. 

Nuclei can be red because of local dust, large foreground 
dust lanes, or because of red stellar populations. One usually 
can distinguish between dust and stellar population redden- 
ing by looking for patchy structures, since they are generally 
due to dust extinction, or to star formation, in combination 
with extinction. In Table 2 we indicate the nature of the 
features that we find in the nuclei. 



The images were calibrated to the STMAG system us- 
ing Holtzmann et al. (1995), using a constant shift of 0.10 
mag to correct to infinite aperture. The method described 
in the same paper was used to iteratively account for the 
colour term in the F450W filter. The colours were corrected 
for galactic extinction using the new dust maps of Schlegel et 
al. (1998), and the Galactic extinction law (Rieke & Lebof- 
sky 1985). After this the bands were K-corrected. This small 
correction (from Persson et al. 1979) was AB — -5z, Al = 
-z and AH = 0. 

To account for the difference in Point Spread Functions 
(PSF) when determining the X — Y colour maps and pro- 
files, the X-band image was convolved with the K-band PSF, 
and the F-band image with the X-band PSF. These PSFs 
were determined with the Tiny Tim package (Krist 1992). 
We then determined minor axis profiles by averaging az- 
imuthally in wedges of 22.5°, centered on the //-band nu- 
cleus. We then combined the surface brightness profiles with 
the ground-based profiles to cover the whole spatial range 
of the bulge, by shifting the ground-based profile on top of 
the HST- profile between about 3" and 6". This way the fi- 
nal profiles have high signal-to-noise everywhere, and the 
high-quality photometric accuracy of HST data, for which 
absolute and relative accuracies of a few percent can be ex- 
pected (Colina et al. 1998). We then subtracted these wedge 
profiles in pairs of bands to obtain colour profiles. In Table 
1 we list the galaxies and their colours at the center and at 
one Reff- 



We find (Table 2) only three galaxies: NGC 6010, 6504 
and 7457, without nuclear dust features, although in some 
cases a large foreground dust lane makes nuclear dust dif- 
ficult to detect. For these three galaxies we analysed the 
isophote shapes more carefully. Dust patches make galaxy 
isophotes irregular in particular the third order Fourier 
terms C3 and S3 (Carter 1978) will be non-zero, and gen- 
erally increasing in amplitude to the blue. S3 and C3 were 
found to be significantly non-zero in both B and / in NGC 
6010 and 6504 but not in NGC 7457. We conclude that the 
signature of dust is to be found in almost all galaxies in the 
sample, including some galaxies indicated with a (-) sign in 
column 3 of Table 2. 

How much reddening is caused by the dust? In Fig. 3 we 
show the colours at the center (filled dots) and at one bulge 
effective radius (open dots) in a / — H vs. B — I colour- 
colour diagram. Fig. 3 (b) and (c) show that in all cases 
the galaxy is redder in the centre than at r e , by sometimes 
very large amounts. Since the vector indicating reddening by 
dust is almost parallel to the the vector indicating changes 
in metallicity, it is not possible to say exactly how much of 
the reddening is due to extinction. Lower and upper limits 
to the internal extinction may be estimated as follows. 

A plausible upper limit to the reddening is derived by 
assuming that the galaxy at 1 R c ff is dustfree. This is per- 
haps reasonable, since at these radii no structure is seen 
in the colour maps, and there is a small dispersion in the 
colour-colour plot. We infer that the central I — H and the 
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Table 2. Nuclear features 



NGC 


F. 


Nuc. 


Blue 


60mu 


lOOmu 


comment 




disk 


d/p 


Patch 


[Jy] 


[Jy] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




5326 


- 


+ 


- 


- 


- 


LGC 361 


5389 


+ 


+ 


- 


0.33 


1.51 




5422 


+ 


+ 


- 


- 


- 


LGC 373 


5443 


- 


+ 


- 


0.33 


1.48 


LGC 373 


5475 


- 


+ 


+ 


- 


- 




5577 


- 


+ 


MDSF 


0.58 


1.94 


LGC 379 


5587 


- 


+ 


- 


0.23 


0.79 




5689 


+ 


+ 


- 


0.46 


1.27 


LGC 384 


5707 


+ 


+ 


- 


- 


- 


LGC 384 


5719 


+ 


? 


- 


8.05 


17.1 


LGC 386 


5746 


+ 


? 


- 


1.33 


8.88 


LGC 386 


5838 


- 


+ 


+ 


0.74 


1.47 


LGC 392 


5854 


- 


+ 


+ 


- 


- 


LGC 393 


5879 


+ 


? 


MDSF 


0.29 


3.04 


LGC 396 


5965 


+ 


+ 


- 


0.39 


1.84 




6010 


- 


- 


- 


- 


- 




6504 


- 


- 


- 


- 


- 




7331 


+ 


? 


- 


23.1 


81.6 


LGC 459 


7457 


- 


- 


- 


- 


- 




7537 


- 


+ 


MDSF 


- 


- 





Notes to Table 2: A '+' sign in column (2) indicates the pres- 
ence of a large foreground disk, in column (3) the presence of a 
red nuclear dust patch or lane, and in column (4) the presence 
of a blue patch near the nucleus. IRAS fluxes are presented in 
column (5) and (6). 'MDSF' in column (4) means 'mixed dust 
and star formation' in the whole nuclear area. The LGC number 
is the group designation (by Garcia 1993). 



B — I colours are reddened by an average of 0.42 ± 0.06 
mag, and 0.61 ± 0.11 mag respectively. A lower limit to 
the reddening can be found if one assumes that the stellar 
populations are never redder than those of the central Virgo 
galaxy NGC 4472. Using the data at 5" (or ~ 400 pc) from 
Peletier et al. (1990) for NGC 4472, converted to B - I and 
I — H as described in the Appendix, we find that the I — H 
colours would be reddened by on the average 0.37 ± 0.06 
mag and B — I by 0.41 ± 0.10 mag. Using the Galactic 
extinction law (Rieke & Lebofsky 1985) these numbers cor- 
respond to an average internal extinction A^ between 0.89 
and 1.01 mag, if I — H is used, and between 0.56 and 0.78 
mag when one uses B — I. The estimates using I — H are 
somewhat higher than those from B — I, because nuclear 
star formation in some galaxies makes the radial colour dif- 
ference in B — I smaller, as is seen in NGC 5838, 5854 and 
5475. In the B — I colour map of these galaxies there are 
blue nuclear features which do not appear in / — H. These 
patches, presumably caused by young stars, reduce the ra- 
dial colour variations in blue colours like B — I but not in 
I -H. 

Our upper limits for the amount of dust inferred here are 
probably somewhat too high, because part of the reddening 
is caused by stellar population gradients. From the most 
dustfree galaxies we infer that the colour difference due to 
stellar populations between center and 1 R ef f is « 0.1 mag 
in both B — I and I — H, which implies that colour gradients 
in the inner parts of bulges are affected much more by dust 
than by stellar populations. 



3.2 Stellar Populations 

Apart from providing high resolution, HST also has the ad- 
vantage that the photometric conditions are very stable, so 
that accurate colours can be determined. For this reason, 
and because we have a very large colour-baseline, we can 
use the colour-colour diagram to infer information about the 
age-spread in the sample, and about the cause of the stel- 
lar population gradients. Studying Fig. 1 carefully, we see 
that many physical phenomena are playing a role in these 
galaxies. We observe the combined effects of extinction, re- 
cent star formation, and old stellar population gradients. To 
disentangle them we have displayed the galaxies in Fig. 1 
according to their morphological type. Although the type 
(from do Vaucouleurs et al. 1991) has been given based 
on low-resolution optical observations, we can see that the 
properties of the galaxies change smoothly as a function of 
type. SO galaxies have small dustlanes, and have rather fea- 
tureless colour maps. Sb galaxies tend to have strong dust- 
lanes, while the Sbc galaxies are considerably bluer, have 
lower surface brightness, show patchy dust and star forma- 
tion together, and are rather different from the rest of the 
galaxies. 

Having established (in Section 3.1) that extinction at 1 
R e ff is probably negligible, we will now consider what the 
colours can tell us about the stellar populations. First we 
note the small scatter amongst the open symbols in Fig. 3a 
confirming that extinction is not important at 1 R c ff. If we 
exclude the three galaxies with the latest Hubble type (the 
open crossed symbols in Fig. 3c) the stellar populations at 
R c ff form a rather tight sequence in the B — I vs. I — H plane. 
Fig. le shows that the central regions of our 3 Sbc are full 
of dust and regions of recent star formation, we will not 
consider them further in this analysis of stellar populations. 
In Fig. 3a we also show a set of Single Stellar Population 
(SSP) models of Vazdekis et al. (1996). The models dis- 
played here have a Salpeter IMF, with a reduced number of 
low-mass stars, to match better the Scalo (1986) IMF. We 
see that, independent of the amount of extinction, the tight- 
ness of the colour-colour relation shows that the luminosity 
weighted age of the stars is very similar from bulge to bulge. 
According to the models the age-spread would be about 1- 
2 Gyr, although the bluest (and generally faintest) bulges 
would be somewhat younger. We derive a mean luminosity 
weighted age of 9 Gyr from these models. While the abso- 
lute ages are poorly constrained, the relative ages are much 
more robust so we can conclude that the age spread amongst 
bulges in this sample is small. (If we use Worthey (1994)'s 
models the inferred luminosity weighted age would be im- 
plausibly small : 2 Gyr). The galaxy sequence runs parallel 
to line of constant age suggesting that the colour-variations 
from galaxy to galaxy are due predominantly to changes in 
metallicity supporting the view that the colour-magnitude 
relation for early-type galaxies (e.g. Bower et al. 1992) is 
mainly driven by changes in metallicity. 

Figure 3c shows that the type dependence along the 
colour-colour relation is very small. The outlying point is 
NGC 5854 which has a blue central patch (probably of re- 
cent star formation). NGC 7457 is a faint SO galaxy that has 
bluer B — R and U — R colours than an elliptical of compara- 
ble luminosity (Balcells & Peletier 1994) accounting for its 
rather young inferred age. The bulk of the galaxies of type 
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Sb and earlier occupy a narrow band in B — I vs. I — H . 
As a comparison we have plotted the colour-colour relation 
for Coma, determined from the data of Bower etal (1992), 
converted to our colours as described in the Appendix. We 
find that the colours of our bulges at 1 R e e are very similar 
to those of bright Coma galaxies. Our bulges however are 
somewhat bluer in I — H and redder in B — /, which, ac- 
cording to the stellar population models, could be explained 
if they are slightly older and somewhat less metal rich. The 
colour-colour conversion is rather uncertain, however. Also 
plotted in Fig. 3c are the colour-colour diagrams converted 
using the theoretical models of Vazdekis et al. (1996) and 
Worthcy (1994). As can be seen in the Figure, the differ- 
ence between the three colour-colour diagrams in B — I and 
/ — K is quite large. Although we argue (see Appendix) that 
our empirical calibration can be trusted much more than 
the theoretical calibrations, one should be careful not to 
over-interpret the data. What seems safe to conclude is that 
Coma galaxies are to be found in the shaded area, and that 
the ages of our bulges are similar to those of the early-type 
galaxies in Coma. 



4 DISCUSSION 

In the previous section we have shown that: 

(i) Centers of bulges of early-type spirals are generally 
dusty. We find that Ay on the average lies between 0.6 and 
1.0 mag, which implies that Ah should be between 0.1 and 
0.2 mag. 

(ii) At 1 R ff Galactic bulges show a very tight I — H vs. 
B — I relation, implying that the age spread among bulges 
of early- type spirals is small (at most 2 Gyr). 

(iii) The colours of bulges of early-type spirals at 1 R c g 
are similar to the colours of early-type galaxies in the Coma 
cluster. This implies that the age difference between nearby 
bulges and cluster-ellipticals is small, probably smaller than 
about 2 Gyr. 

Since the extinction observed, is more than a factor 2 
(0.76 mag), if the dust is located in or near the center we can- 
not see through to the other side of the galaxy. The extinc- 
tion locally will be very large (see Sadler & Gerhard 1985). 
The fact that this occurs in so many of our galaxies shows 
that the very central region is almost always optically thick 
in B or F. The situation with dust in the centers of spiral 
bulges is similar to that of elliptical galaxies. Van Dokkum 
& Franx, analyzing WFPC data of nearby elliptical galaxies, 
found evidence for central extinction in 75% of their sam- 
ple. This is more than was previously found from the ground 
(Ebneter et al. 1988, Veron-Cetty & Veron 1988). Goud- 
frooij et al. (1994b) detected dust in 41% of the galaxies 
in their sample of Revised Shapley-Ames (Sandage & Tam- 
mann 1981) galaxies. A similar detection rate as found from 
IRAS fluxes by Knapp et al. (1989). The higher angular 
resolution of HST allows us to reach much lower detection 
limits, we have been able to detect dust in 95 ± 5% of our 
sample of bulges of early-type spirals. This detection rate is 
higher than for ellipticals, although our method of finding 
dust, using colour images, is more sensitive than the method 
of van Dokkum & Franx, who determined their dust masses 
from one V-band image only. Van Dokkum & Franx found 



an average dust mass of 4 x 10 3 Mq . For the bulges analysed 
here we find an average dust mass of ~ 10 4 Mq, determined 
using the method of van Dokkum & Franx: for each dust fea- 
ture the mass is calculated using Md us t = £ (Ay) T^ 1 , with 
E the area of the feature, (Ay) the mean absorption in the 
area, and Ty the visual mass absorption coefficient (Sadler 
& Gerhard 1985). Ty is taken to be 6 10~ 6 mag kpc 2 M" 1 . 
Rough values of (Ay) were determined by taking the differ- 
ence in B — I or I — H of the feature and the dustfree values 
at 1 R c ff , converted to Ay using the Galactic extinction law 
(Rieke & Lebofsky 1985). If we assume a Galactic gas to 
dust ratio of 130 (see van Dokkum & Franx 1995) we find 
the galaxies analysed in this paper have an average of 10 6 
Mq of interstellar material in their nuclear regions. 

The origin of the nuclear dust is unclear. At larger radii 
kinematic observations of gas, usually associated with dust, 
in ellipticals show that it is often decoupled from the stellar 
velocity field (for a discussion see Goudfrooij et al. 1994b). 
This is used to imply that large-scale dust lanes are of ex- 
ternal origin perhaps being accreted during galaxy mergers 
or interactions. However, the origin of the small arcsecond- 
scale dust lanes found in ellipticals, often oriented along the 
major axis of the stellar body (Goudfrooij et al. 1994b), 
could well be internal. Scaling from the numbers for stellar 
mass-loss for bright ellipticals given by Faber & Gallagher, 
1976, we estimate that for typical bulges in our sample the 
mass will be dposited into the bulge at the rate of « 0.1 - 1 
Mq yr _1 , so there is no problem accruing the dust we see. 
The amount of dust, the large detected fraction, and the 
fact that the dust lanes are found parallel to the major axis 
indicate an internal origin. The central dust provides a suit- 
able environment for centrally concentrated star formation, 
probably leading to strong metallicity enhancements in the 
central 100 pc of bulges. In that regime, it is possible that so 
much dust is produced that ordinary dust destruction mech- 
anisms are ineffective. Potentially, high resolution spectral 
observations will show us whether this dust is indeed of in- 
ternal origin, and what metallicities that are being reached. 
The dust however could provide some serious observational 
difficulties for the determination of the inner slope of the 
stellar surface density profile of bulges, even in the H band. 

Several of our bulges have central features which resem- 
ble the inner disk in the giant elliptical NGC 4261 (Jaffe et 
al. 1996). Examples are NGC 5326, NGC 5587, NGC 5838 
and NGC 5854. NGC 4261 has a LINER spectrum in the 
center, Jaffe et al. argue that this inner disk might provide 
the fuel for the AGN. Only 4 of the objects studied in our 
paper however are known to be (weakly) active, NGC 5746, 
NGC 5838, NGC 5879 and NGC 7331, all classified by Ho 
et al. (1997) as T2, transition objects between LINER and 
Seyfert, with narrow emission lines. Their colours are how- 
ever entirely consistent with an extincted stellar population. 
Furthermore, many of the ellipticals with similar features are 
not active galaxies. It seems that an inner disk of gas and 
dust alone is not sufficient to produce an AGN. 

No HST colour profiles of bulges of spirals have ap- 
peared in the literature up to now. The only paper present- 
ing high-resolution HST colour profiles of ellipticals is Car- 
olio et al. (1997a), who report V and I profiles of 15 galax- 
ies with dynamically decoupled cores using the refurbished 
WFPC2. They find that their galaxies all have very similar 
V — I gradients between radii of 1.5" and 10", while the dis- 
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persion is larger between 0.25" and 1.5", where some galax- 
ies are seen with gradients that are twice as large, while some 
others are much smaller, or even negative. Clearly the be- 
havior of V — I is different in their inner 1.5" as compared to 
the area further out. Although Carollo et al. (1997a) masked 
out dusty areas before obtaining the colour profiles by radi- 
ally averaging the remaining light, this process will proba- 
bly not have removed all the extinction. The large spread in 
colour gradient in the inner region might be a confirmation 
of significant quantities of dust near the center in almost all 
elliptical galaxies. 

What can we learn about the formation of bulges? The 
fact that the ages of most of the bulges in this paper are so 
similar and old makes it very difficult to form the bulges of 
early-type spirals (SO-Sb) through secular evolution of disks. 
In this scenario it is expected that bulges regularly undergo 
major bursts of star formation, to convert gas that has been 
tunneled to the central area through the presence of a bar 
into stars. This would mean that we would expect to find 
more young bulges and a large spread in bulge ages. Another 
problem for early-type bulges is that the bulge densities are 
up to a factor 100 higher than in the center of the disks. 
To create such high overdensities in the disk, the required 
star formation rates are such that they could easily disrupt 
the disk (Ostriker 1990). Late type galaxies (Sbc's and later 
types) might be different. We see that the three in our sam- 
ple have younger ages (crossed circles in Fig. 3), and also 
the density contrast between bulge and disk is much smaller. 
For these galaxies in the bulge region spiral arms, blue star 
formation regions and dust lanes are seen, as in the disk, 
telling us that the stellar populations of bulge and disk here 
are very much the same (see Fig. le). Although Norman 
et al. (1996) use the fact that the Fundamental Plane has 
changed little from a redshift of the order of 0.5 to now as an 
argument against secular-evolution driven bulge building, it 
is not clear whether bulges of late-type galaxies really lie on 
the fundamental plane of ellipticals and early-type bulges, 
since no data is available at present. 

The fact that the ages of our bulges are all so simi- 
lar supports the idea of a major episode of star formation 
in the past, in which most of these bulges (and also the 
bright galaxies in large clusters like Coma) were formed. It is 
currently thought that cluster ellipticals must have formed 
at redshifts beyond z — 3 (Ellis et al. 1997, Stanford et 
al. 1997), because of the lack of evolution seen in clusters 
of intermediate redshift. This corresponds to an age of 10.5 
Gyr (fi = 0.2, 0\ = 0.8, H = 70 km s" 1 Mpc -1 ), or 
8.4 Gyr (fi = 1, fi A = 0, H„ = 70 km s" 1 Mpc" 1 ) (Hogg 
1999). The fact that the colours of the majority of our bulges 
are similar to the Coma galaxies indicates that our bulges 
(except maybe the Sbc's) are also old, and formed at red- 
shifts beyond z — 3. These observations endorse the original 
model by Eggen et al. (1962) forming bulges in the begin- 
ning during a monolithic, or clumpy, collapse. The fact that 
the colours of bulges and inner disks are very similar (Tern- 
drup et al. 1994, Peletier & Balcells 1996) then implies that 
the disk formed very gradually from inside out, with the 
age of the inner disk similar to the age of the bulge. Can 
we then say something about the age of the disk? Would 
it be possible that the whole disk formed at the same time 
as the bulge? At that point a conflict arises with the star 
formation history of the universe derived e.g. by Madau et 



al. (1996) on the basis of data of the Hubble Deep Field, 
which shows a maximum between z — 1 and 3. Since our 
bulges are found in a variety of environments, they are in no 
way special, and if most early-type spirals would form at red- 
shifts beyond 3 this would also imply that the maximum in 
the HDF would have to go to larger z. However, this prob- 
lem would be solved if disks of early-type galaxies on the 
average would be considerably younger than their bulges. 
With a sufficiently large age difference between bulge and 
disk the luminosity-weighted age of the galaxy (which is in 
general dominated by the disk) can then be young enough. 
There is nothing in the data of Terndrup et al. 1994 and 
Peletier & Balcells 1996 preventing this from happening. 

Alternatively, Kauffmann (1995) and Governato et 
al. (1998) have pointed out that in biased models of galaxy 
formation accelerated evolution can be expected in dense 
regions. Our observations here show that there is no differ- 
ence in age between the ellipticals in Coma and our galactic 
bulges. This would mean that our bulges, except for the lat- 
est types, must have formed early, beyond z — 3, as well. On 
the other hand, we also do not find any dependence of bulge 
age upon environment. The bulges are found in all kinds of 
environments, ranging from isolated to groups of more than 
20 members (see Table 1). Apparently everywhere in the 
universe these intermediate-size galaxies must have started 
forming early-on. 

It appears that the only feasible solution is that bulges 
of early-type spirals are old, and disks considerably younger. 
This is in agreement with Abraham et. al. (1998), who find 
for the disk galaxies in the HDF that their bulges are signif- 
icantly older than their disks (up to 50%). One has however 
to take into account that at z ~ 1 many large spiral galaxies 
are found (Lilly et al. (1998), whose stellar populations are 
consistent with a declining activity since at least z=1.5 - 
2 (Abraham et al. 1998). Semi-analytic galaxy formation 
models also support this picture. In Fig. 4 a histogram is 
shown of the distribution of V-band averaged ages of ellip- 
ticals and bulges in the simulations of Baugh et al. (1996) 
(Baugh, private communication). The models find that the 
large majority of bulges is old. This is independent of the 
environment, contrary to elliptical galaxies, which are found 
to be old in rich clusters, and can be 5 Gyrs younger in the 
field. Bulges are old because the accompanying disk needs 
time to form without being disrupted. 

We have found that bulges of early-type spirals in gen- 
eral are old, that their age-spread is smaller than 2 Gyr, 
and that their colour gradients are mainly due to metallicity 
gradients in the stellar populations. These results once more 
support the idea that bulges and ellipticals are similar ob- 
jects. The good agreement between the observational results 
of this paper and models like the semi-analytical models of 
Baugh et al. (1996) strongly supports the picture that bulges 
were formed through monolithic or clumpy collapse, and 
formed much before theier disks. However, one can see that 
bulges of Sbc galaxies are different in many respects, and it 
is likely that these differences will be larger for Sc galaxies 
and later types (see also the recent HST-study of Carollo 
et al. (1997b, 1998). They are smaller, younger, have lower 
central surface brightness, and it is not clear whether they 
fall on the fundamental plane of elliptical galaxies. Since 
the observations until now are not good enough to establish 
whether these late type bulges can be made from bars, and 
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in this way the galaxy type can be changed by secular evolu- 
tion, it is clear that it is very important to study the bulges 
of late type spirals to understand the formation of bulges. 
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Figure 1. B — I and I — H colour maps, and intensity-maps of the galaxies. Greyscale levels for B — I are 1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.7,3.0 
and 3.5 mag, and for I — H 1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4 and 3.0 mag. Superimposed on the colour maps are H-band contours with 
levels from 20 to 10 mag with intervals of 1 mag. The B — I — H maps arc black and white representations of 'real colour' maps made 
from B, I and H, and arc uncalibratcd. a: SO galaxies. 



Figure 1. b: SO/a galaxies. 



Figure 1. c: SO/a-Sab galaxies. 



Figure 1. d: Sab-Sb galaxies. 



Figure 1. e: Sb-Sbc galaxies. 
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Figure 2. Individual colour profiles for the sample galaxies. Each galaxy is represented by a line from the center to the point where 
the bulge contribution is as large as that of the disk. Morphological types are indicated by line types: T < - drawn lines; T = 1: long 
gashed; T = 2: dot - long dashed; T = 3: dot - short dashed; T = 4: dotted. 
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Figure 3. Colour-colour diagrams for the 20 galaxies. In Fig. 3(a) and (b) are displayed the positions of the center (filled) and 1 
bulge effective radius on the minor axis (open circles). Reddenings vector for a reddening of Ay = 0.5 mag (3b) and 1 mag (inset of 
3a) are given as well. Superimposed in Fig. 3a are SSP models by Vazdckis et al. (1996). Solid lines are lines of constant metallicity, 
dashed-dotted lines are loci of constant age. In Fig. 3c the same galaxies at 1 R e ff are plotted, with their symbols coded as a function 
of morphological type. Also added are NGC 4472. An attempt has been made to convert the data of Bower et al. (1992) to B — I and 
I — H, using the models of Worthcy (1994, long-dashed), Vazdckis et al (short-dashed) and the empirical calibration described in the 
Appendix (drawn line)). Since these calibrations do not agree very well with each other we can only say that the early-type galaxies of 
Coma are to be found in the shaded area, and that the position of our bulges is consistent with their colours. 
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Figure 4. Mean V-band luminosity weighted ages of ellipticals and bulges in clusters and in the field, predicted by semi-analytic 
simulations by Baugh (private communication), following prescriptions given in Baugh et al. (1996). Drawn lines are bulges, while 
elliptical galaxies are indicated by dashed lines. 
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Figure 5. V — I - My colour-magnitude relation for the early-type galaxies of Goudfrooij ct al. (1994a) (astcriscs) together with the 
data of Bower et al. (1992) converted to V — I (data for Coma are shown by filled squares and for Virgo by open squares). See text for 
more details 
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APPENDIX A: COLOUR TRANSFORMATIONS FOR THE COLOUR-COLOUR DIAGRAM OF THE 
COMA CLUSTER 

In this appendix we describe how a colour-colour relation for the Coma cluster was obtained in B — I and I — H. We started 
with the excellent data of Bower et al. (1992) in U, V and K. As a first attempt we tried to make the colour conversions 
using single burst theoretical models for old stellar populations using the Salpeter IMF. A fit to the models of Vazdekis et al. 
(1996) gives 

I-K = 0.546 (V - K) + 0.275 
and 

B-V = 0.433 (U-V) + 0.296 

I — K was finally converted to I — H using the new GISSEL96 models of Bruzual & Chariot (see Leitherer et al. 1996) 
using 

I-H = 0.889 (I-K) + 0.034 
since the Vazdekis' models do not tabulate the //-band. Alternatively, a fit to Worthey (1994)'s models gives: 

I-K = 0.707 (V-K) - 0.253 
and 

B-V = 0.407 (U-V) + 0.379 
using a very similar transformation from I — K to I — H: 

I-H = 0.895 (I-K) + 0.047 

There are several reasons to be wary of this approach. In the first place one knows (e.g. Chariot et al. 1996) that the 
infrared colours of the Worthey models for large metallicities are too red by as much as a magnitude. Secondly, the model 
colours of Vazdekis et al. in the infrared, just like those of Bruzual & Chariot (Leitherer et al. 1996) do not appear to be very 
accurate either, since they show abrupt jumps as a function of age and metallicity, due to the rather unknown contribution 
from stars in their late stages of evolution. The conversion between I — K and I — H is only of minor importance. When one 
compares the results of both methods, one finds that the difference in the colour-colour diagram of Coma, when using two 
independent models, is quite large (see Fig. 3c). This shows that a different, independent way to convert the colours would 
be very useful, and for that reason we attempted to derive an empirical conversion. 

In the search of data in B or / of the Coma cluster to combine with Bower et al.'s dataset, we have found the data of 
J0rgensen et al. (1994), with photometry in V and Gunn r, and the Ph.D. thesis of Steel (1998), with V and R. Since there is 
no /-band photometry available, we decided to determine a conversion from V — K to V — R, a band with information very 
similar to V — I, and then convert V — R to V — I using published stars. It was also found that if one determines a colour of 
a galaxy using its integrated magnitudes, the observational error would be much larger than if the same aperture was used 
both times. For that reason we preferred to use the V — R data inside 20" of Steel (1998), for which the dispersion in the 
V — R vs. Rt was only 0.02 mag, rather than the total colours of j0rgensen etal. (1994). Combining the data of Steel (1998) 
with those of Bower et al. (1992) one obtains 

V-R = 0.336 (V-K) - 0.501 

with a scatter in V — R of 0.03 mag. To convert V — R to V — I we fitted a least squares relation to all the Landolt (1992) 
standard stars. We find a very tight relation 

V -I = 2.029 (V-R) - 0.018 

with negligible scatter. The validity of this conversion between V — K and V — I can be for example be established 
by looking at the V — I - magnitude relation for the elliptical galaxies of Goudfrooij et al. (1994a). This comparison gives 
satisfactory results. Finally, we need to obtain a relation between U — V and B — V . Since here also there is quite a difference 
between the theoretical models, we have determined a least squares relation between the observed U — V and j0rgensen's 
B — r, which we first converted to B — V using the conversion given in j0rgensen (1995): 

© 0000 RAS, MNRAS 000, 000-000 



14 Peletier et al. 

B-V = 0.673 (B-r) + 0.184 
to get: 

B-V = 0.356 (U-V) + 0.448 

Errors are difficult to determine. We can get an estimate by comparing the data of Goudfrooij et al. (1994a) with the 
V — K data of Bower et al. converted to V — I (Figure 5). Plotted on the x-axis for the data of Bower et al. are absolute 
F-band magnitudes, obtained using a distance modulus of 30.82 for Virgo and 34.51 for Coma (Aaronson et al. 1986), to 
which a correction of 0.76 mag has been applied to convert them to integrated magnitudes. For the data of Goudfrooij et al. 
(1994a) we plot their My values. There are three galaxies in common between the 2 samples, for which the My values agree 
within 0.10 mag. Comparing the two samples we find that for the galaxies between My = -21 and -24 V — I is redder by 
0.032 mag on the average in the data of Goudfrooij et al. (1994a). This shows that our conversion from V — K to V — I is 
probably reasonable, and based on this, and on the tightness of individual relations, we believe that the error in V — I is at 
most 0.06 mag, the same for I — H or I — K, and 0.05 mag in B — V , resulting in 0.08 mag in B — I. We not only applied 
these transformation to the mean relation for the Coma cluster, but also for the giant elliptical NGC 4472 (see Fig. 3). 



© 0000 RAS, MNRAS 000, 000-000 



This figure "figla.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/99 1 1 5 3 v 1 



This figure "figlb.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/99 1 1 5 3 v 1 



This figure "figlc.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/99 1 1 5 3 v 1 



This figure "figld.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/99 1 1 5 3 v 1 



This figure "figle.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/99 1 1 5 3 v 1 



